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Abstract 

We present detailed discussions on a new approach we proposed in a previous 
paper to numerically study quantum spin systems. This method, which we will call 
re-structuring method hereafter, is based on rearrangement of intermediate states 
in the path integral formulation. We observed our approach brings remarkable im- 
provement in the negative sign problem when applied to one-dimensional quantum 
spin 1/2 system with next-to-nearest neighbor interactions. In this paper we add 
some descriptions on our method and show results from analyses by the exact di- 
agonalization and by the transfer matrix method of the system on a small chain. 
These results also indicate that our method works quite effectively. 
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Section 1 Introduction 

Recent development in experiments on condensed matter has shown us very in- 
teresting systems where strong quantum effects should be realized. One of them is 
the quantum spin system in low dimensions. Much theoretical work has been done 
to obtain quantitative as well as qualitative properties of these quantum effects [|TJ . 

One powerful tool to numerically investigate quantum spin systems is Monte 
Carlo approach using the Suzuki- Trotter formula [@] . Study through this method has 
brought us very intriguing results on the ferromagnetic system. It is widely known, 
however, if one applies this method to frustrated systems one often encounters the 
so-called negative sign (NS) problem [[|. The NS problem, which becomes more 
serious on larger lattices, makes it very difficult to get statistically meaningful results 
in numerical calculations. 

In a previous paper we showed that a new approach (re-structuring (RS) 
method) is useful to obtain Monte Carlo results for one-dimensional quantum spin 
1/2 system with next-to-nearest neighbor interactions, which suffers from the NS 
problem in the antiferromagnetic case. Essential point of RS method is to choose 
a set of states for the path integral which is appropriate to numerical calculation. 
Conventional choice of the set, whose states consist of eigenstates of ^-component of 
the Pauli matrix on each lattice site, is the simplest one but the NS problem turns 
out to be serious with this set. Since any other choice is possible as long as the 
set is complete, we employ a set made of eigenstates of local Hamiltonian for every 
two neighboring sites. Although a complete solution for the NS problem could not 
be obtained, our results, which show much improvement in numerical calculations, 
inspires us with confidence that quantum Monte Carlo method effectively works if 
a set of states for the path integral is appropriately chosen. 

In this paper we add descriptions on our method in some detail. We then analyse 
the same system by the exact diagonalization and by the transfer matrix method 
on a small chain. We will see these analyses are also helpful to confirm superiority 
of RS method to the conventional one. In section 2 the model and the conventional 
approach, which is necessary in our analyses for comparison, are scanned. Section 3 
is for detailed descriptions of RS formulation. Results from our analyses are given 
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in section 4 and final section will be devoted to summary and discussions. 

Section 2 Model and Conventional approach 

The system we study is the quantum spin 1/2 system with next-to- nearest neigh- 
bor interactions on a one-dimensional chain, the simplest one among those suffering 
from serious NS problem. The Hamiltonian of this system is 

1 N 

H = - + 0^+2), (1) 

1 i=i 

where N is number of sites on the chain and &N+i = o% (periodic boundary condi- 
tion). The partition function Z is given by Z = trie' 1311 ) with inverse temperature 
f3. 

Let us describe the conventional approach. State on each site is represented by 
^-component of the spin, namely up and down, or + and — . In this representation 
states of the system are given by 

I OL > = \ S 1 ,S 2 ,S 3 , ...,S N >, 

where Sj = + or — . The identity operator is then 

1 = X] I sx,s 2 ,—,sn >< s 1} s 2} ... } s N I . 

To use this identity operator in the Suzuki- Trotter formula, we divide the Hamil- 
tonian into four partsQ as schematically shown in Fig. 1(a). Thus we come to an 
expression 

Z = lim tr{(e-^ l/n e-^ 2/n e-^ 3/n e-^ 4/ Ti, (2) 

n— >oo 

where 

1 N/2 

Hi = - <?2i-i<?2i, 

1 1=1 

1 N / 2 

H 2 = - ^2i<?2i+l, 

1 1=1 

1 For techincal reasons number of sites in this case should be restricted to be multiple of four. 
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1 N/i 

1 i=l 
1 N/i 

Hi = — y^(0"4i-l<74t_)_l + <?4j<?4i+2)- 
i=l 



With above complete set and partial Hamiltonians we obtain the partition function 
used in Monte Carlo calculations, 

z c ) = w{{a h a' p a",a'"}), (3) 

{a u a'.,a" ,a"'} 

where 

x < «J | e"^/™ | aj' >< af | e"^ 4 ^ | a J+1 >, 

suffix j numbering sites along the Trotter axis and a n +i = «i- 

In this system w({aj, a'j, a", a"'}), total product of expectation values over one 
configuration, can be negative. Appearance of this negative weight brings the NS 
problem. In order to numerically calculate some physical quantity (A) one then 
should subtract contributions of negatively signed configurations, A_, from those of 
positively signed ones, A + . Namely, 

^ = Y^t- « 

where Z + (Z_) is number of configurations with positive (negative) weight. The 
result would suffer from serious cancellation when Z_ ~ Z+. 



Section 3 RS Method 

Formulations of RS method is outlined in our previous paper |4j] . Here we briefly 
repeat them and add detailed information on its effective Hamiltonian. 
We start from rewriting ([]]) in the following form 

1 N/2 

H — — ^(<?a,i<?a,i+l + <?f>,i<?b,i+l + &a,i&b,i + <3b,i& a,i+l ) , 
i=l 



where we denote odd and even sites with suffix a and b, respectively, 

&a,i = 02i-l) ®b,i = 0"2i- 

We employ the complete set for which the operator products (T a ,i&b,i are diagonalized. 
Explicitly, we use 

| a >=\ Si, S 2 , S N/2 >, (5) 
where Si = lj, ©i, 0j or — lj with 

I lj ^* | ~t~o,t) ~l~6»i -^j 
I ©i > = 7^(1 "t"a,zj — 6,i > ~KI ~~ a,ii +6,i >)? 

I ©i > = 7^(1 ~t~o,i) — 6,i > — (| ~~ o,ij +6,i >)) 
| lj ^* | o,i; 6,i ^* • 

Hamiltonian (|]) is accordingly devided into "odd" and "even" parts 

H = H + H e , (6) 

where, as shown in Fig. 1(b), 

i W 

h =-y: hoi (7) 

z i=i 

hoi — {&a,2i-l&a,2i + &b,2i-l&b,2i + &b,2i-l&a,2i 
+ 2 <T a,2i-l°"6,2j-l + 2 <T a I 2iO"6,2iJ, 

i iV/4 

#e = ^E ^, (8) 
Z i=l 

ft-ei = (<?a,2i<?a,2i+l + &b,2i&b,2i+l + &b,2i&a,2i+l 

Expression of the partition function with these partial Hamiltonians is 

Z = lim tr{(e-^ o/n e-^ e/n ) n }. (9) 



Inserting identity operators made of a's in ([5]) between the exponents we obtain 
partition function in RS method, 

n 

Z RS = E II < a 3 I e ~ Pil °' n I a 'j >< a 'j I I a 3+l > • ( 10 ) 

K, Q ;.}i=i 

In order to obtain Boltzman weights to be used in Monte Carlo simulations, we need 
to construct < S' 2i _ x ,S 2i \ exp(—^h oi ) | S 2i ^i,S 2i > from local matrix elements 
< S^i-iJ^i I h a I S 2 i-i,S 2 i > as well as < S' 2i ,S' 2i+1 \ exp{-^h ei ) \ S 2i ,S 2i+ i > 
from local matrix elements < S' 2i , S' 2i+1 \ h ei \ S 2i , S 2i+ i >. Since these matrix 
elements are ^-independent and 

^ S 2 i-n S' 2 i I hoi I S 2 i-i, S 2 i >=< S 2i , S 2i+1 | h e i \ S 2 {, S 2 j+i >, 

we only need to know values of 

< S[, S' 2 | h \ | Si, S 2 >=< S[, S' 2 | <7 0j i(7 aj 2 | Si, S 2 > 

+ < S[, S 2 | &b,\&b,2 I S'2 > + < Si, S 2 | <Tfe,l<T a ,2 I 5*1, 5*2 > 
+ - < S[, S 2 I (Ja.loVl I Si, S'2 > +2 < Si, S 2 I f7 0)2 (7&,2 I Si, S 2 > . (11) 

It is an easy task to calculate each matrix element in (|TT|). For example, 

< -li, 1 2 I (T a ,l&a,2 I ©1,02 > 
= < -a,l, -6,1 |< +a,2, +6,2 | (2o"Ji(7~ 2 + 2(7~ 1 0-+ 2 + e^jO^) 
7/^(1 +0,1 J ~ 6,1 > + I — o,l j +6,1 > )^/2 = (l +a,2, — 6,2 > — I — a,2, +6,2 >) 

=< — li, 1 2 I (| li, -I2 > - I -li,l 2 > + I 01,02 >) = -1, 

where 

a ± = X -{o x ±io v ). 
Values of < S[, S 2 \ h Q i \ Si, S 2 > are summarized in Table. 
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Section 4 Transfer Matrix and Exact Diagonalization 

In this section we present results from analyses with RS partition function (10) 
on an N = 8 chain by the transfer matrix method H] and the exact diagonalization. 
Results from conventional partition function (3) are also provided for comparison. 
Purpose of these analyses are twofold. One is to confirm RS algorithm and previous 
Monte Carlo results. Another is to study how systems with Zq and Z^g depend 
on Trotter number n. 

First let us show results for ratio of negatively weighted configurations to total 
configurations, P = Z_/(Z + + Z_), calculated by the transfer matrix method. We 
use 'real-space' transfer matrix to realize the same boundary condition as those in 
other methods. In the conventional approach we calculate product of matrix Tk, 
whose matrix element (Tk) a 'a corresponds to < a' | exp(—/3Hk/n) | a >. Therefore 

Z£ ] =tr((T 1 T 2 T 3 T,) n ). 

Similarly with (T ^) a > a =< a' \ exp(—(3H ( e )/n) \ a >, 

Z$=tr((T e T o r). 

In figure 2(a) we plot P from the transfer matrix method together with previous 
Monte Carlo results. We see Monte Carlo data and transfer matrix results for n = 2 
agree very well. Also we observe agreement between those for larger Trotter numbers 
is satisfactory^. Our Monte Carlo work is thus nicely supported by an analytical 
calculation. 

Next we calculate 

= Z + - Z_ = Z_ 

Z + + Z_ Z'' 1 ; 

which equals to 1 — IP and is the ratio of the partition function Z obtained from the 

Hamiltonian ([!]) to the partition function Z' calculated from the modified Hamilto- 

nian that gives absolute weights 0. 

For iV = 8 chain it is not difficult to calculate all energy eigen values Ei (i = 

1, 2^) by the exact diagonalization. We therefore can obtain value of exact par- 

2 Small discrepancy in low temperature region might suggest that additional improvements in global flips, which 
are related to the ergodicitylfj , are necessary. 



tition function at any temperature by equation 



(13) 



The same is true for the denominator in fll2|) ; 



Z' = J2exp(-(3E[), 



(14) 



i=l 



with energy eigen values E\ for the modified Hamiltonian. We refer to the ratio 
R calculated from (p~3f ) and ( |14"D as -Rdg- It should be noted that value of Rdg 
depends on method used to calculate the denominator Z'. We will see RS method 
gives better ratio compared to the conventional method. 

Before going to full numerical calculations, let us comment on ground state en- 
ergy of each Hamiltonian. As discussed in ref.[Q, at very low temperature where 
contribution from the lowest energy eigen value is dominant, 



because Z + — Z_ ~ exp(—(3Eo) and Z + + Z_ ~ exp(—/3E' Q ) with E (E' Q ) denoting 
energy of ground state for the Hamiltonian (for the modified Hamiltonian). Exis- 
tence of severe NS problem implies E > E' . 

Our results are in agreement with the above discussion. Values we obtain by 
the exact diagonalization are E = —8.25 and E' Q = —10.14 (—11.32) for RS (con- 
vensional) method. The fact that the energy difference Eq — E' in RS method is 
smaller than that in the conventional method indicates improvement introduced by 
RS method. 

For the second purpose we also calculate R by the transfer matrix method with 
several values of Trotter number n. All numerical results on R are summarized in 
figure 2(b). Here we would like to point out two features. One is that, as is expected, 
the Rdg is much larger in the RS method than in the convensional method, especially 
for low temperatures. Another feature is rather unexpected one. In RS method it 
turned out that i?'s for finite Trotter numbers, which should coincide with Rdg i n 
n — > oo limit, are always larger than Rdg- 111 the conventional method, on the 
contrary, all R values calculated by the transfer matrix method lie below Rdg- This 



R ~ exp(-(3(E - E' Q )) 
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encourages us to promote RS method in Suzuki- Trotter formula because it becomes 
possible to obtain statistically meaningful Monte Carlo results even for small values 
of n. 

Finally we present results for the energy of the system. In figure 3, exact values 
obtained by the exact diagonalization, 

nJV 2 N 

E = J2E l .e-^/J2e^ E % 
i=i i=i 

and other values calculated by the transfer matrix method are plotted. It is clear that 
conventional method fails to give reliable values when Trotter number is small. Even 
when n = 8, difference from exact one becomes large as temperature T decreases. In 
the RS method, on the contrary, we get satisfactorily stable results even for n = 2 
and all data seem to nicely converge to the exact ones when n is enlarged. 

Section 5 Discussions 

In this paper we provide additional explanations of RS approach, which we pro- 
posed in a previous paper j|] in order to make meaningful Monte Carlo study of 
quantum spin systems suffering from the NS problem. From technical point of view 
main task is to analytically calculate matrix elements for effective Hamiltonian in 
( |lOl) . Descriptions on it are presented in section 3. We also present results from anal- 
yses of quantum spin 1/2 system with next-to-nearest neighbor antiferromagnetic 
interactions on a small chain, applying RS approach in the transfer matrix method 
and the exact diagonalization. In section 4 we see all results obtained through these 
analyses indicate that RS method is effectively works in this case. 

Let us comment on possible applications of RS method to other quantum spin 
systems. Qualitatively the RS method will be effective for any frustrated systems 
suffering from the NS problem. But its power to bring how much quantitative 
improvement in numerical calculations does strongly depend on system's interactions 
and spacial constructions of clusters. 

One interesting application would be to investigate quantum spin 1/2 XY chain 
with next-to-nearest neighbor interactions. For qualitaitve discussion we calculate 
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energy eigen values of the system on iV = 8 chain by the exact diagonalization. The 
lowest enegy is Eq = —5.78 for the XY Hamiltonian while E' = — 7.14(— 9.44) for 
the modified Hamiltonian in RS (conventional) method. These results suggest RS 
method works more powerfully on the XY chain in comparison with the Heisenberg 
chain, which should be one example of the above statement. 

Another intriguing application was made to the A chain by Nakamura and Saika. 
They reported that RS method is quite useful for this system 0. 

We would like to point out a very lucky case where the NS problem is completely 
solved, although it might be less interesting from physical point of view. This is 
one-dimensional antiferromagnetic quantum spin 1/2 system with next-to-next-to- 
nearest neighbor interactions in addition to next-to-nearest neighbor one, where all 
coulpings have the same strength^. This system sufferes from severe NS problem in 
the conventional approach. When one chooses eigen states of the nearest interaction 
for the nearest two spins, matrix elements of the effective Hamiltonian become quite 
simple so that negatively weighted configurations do not appear at all. 

In two-dimensional cases systems on triangular lattice will be most important in 
application. If we employ eigen states for cluster of two neighboring spins as we did 
in one-dimensional case it is not unique which two should form a pair because more 
than one neighboring sites exist on the lattice. Symmetry of the lattice inspires us 
clustering of three spins is appropriate, for which one has to do with more compli- 
cated effective matrix elements. At present it is an open question which way gives 
better results. 

In any case a good choice of complete set which improves the NS problem will 
help us to gain deep insight into the ground state of the system under investigation. 

3 Number of sites of this system also should be multiple of four for technical reasons. 
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Table Caption 

Values of matrix elements < S[,S' 2 \ h a \ \ Si, S 2 > described in section 3. Rows 
and colomns denote states < S[, S' 2 | and | Si, S 2 >, respectively. All other matrix 
elements are zero. 



Figure Captions 
Figure 1 

(a) Schema to show how to divide the Hamiltonian ([!]) into four parts on an eight- 
site (N = 8) chain with periodic boundary condition. Sites are numbered from one 
to eight. Filled circles, squares, diamonds and triangles on links denote interactions 
to belong partial Hamiltonians Hi, H 2 , H3 and in (|2]), respectively. 

(b) Schema to show how to re-divide the same Hamiltonian ([[]) into H a and H e 
in (H). Sites are renumbered from al to 64. Open diamonds (triangles) denote 
interactions to fully belong to partial Hamiltonian H Q (H e ). Interactions on links 
with open circles should be equally shared by H Q and H e . 

Figure 2 

(a) Ratio of negative weight calculated on an iV = 8 chain by Monte Carlo simu- 
lations (MC) and by the transfer matrix method (TM), using the RS Hamiltonian 
(RS) or the conventional one (conv.), with several values of Trotter number n. 

(b) Ratio of partition functions on an iV = 8 chain. The exact diagonalization (DG) 
is employed in addition to the transfer matrix method. 

Figure 3 

Energy per site calculated on an N = 8 chain by means of the exact diagonalization 
(DG) and the transfer matrix method (TM). The transfer matrix method is applied 
to both Hamiltonians (RS, conv.) with Trotter number n =2, 4 and 8. 
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